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We compute analytically the dominant contribution to the tree-level bispectrum in the Starobin- 
sky model of inflation. In this model, the potential is vacuum energy dominated but contains a 
subdominant linear term which changes the slope abruptly at a point. We show that on large scales 
compared with the transition scale fco and in the equilateral limit the analogue of the non-linearity 
parameter scales as (k/kof' , that is its amplitude decays for larger and larger scales until it becomes 
^ ' subdominant with respect to the usual slow-roll suppressed corrections. On small scales we show 

' that the non-linearity parameter oscillates with angular frequency given by 3/fco and its amplitude 

, grows linearly towards smaller scales and can be large depending on the model parameters. We also 

compare our results with previous results in the literature. 

^: 

^ • I. INTRODUCTION 

as : 

' The idea of inflation in the very early universe [iHll is a very compelling solution to the problems of the Big Bang 
__, theory. However one should keep in mind that there may be other viable alternative explanations. Other successful 
' predictions [p of the simplest inflationary models are the nearly scale invariant, mostly of adiabatic origin and nearly 
r \ ' Gaussian primordial curvature perturbation. This primordial perturbation is generated due to quantum vacuum 
fluctuations of the scalar field driving inflation and will at a much later stage in the evolution of the universe give 
origin to the cosmic microwave background (CMB) radiation anisotropics and the large-scale structure of galaxies in 
the universe. 

Q ' While the simplest models are appealing due to their simplicity and fewer ingredients they are not always realistic 
and often more realistic models introduce several fields, non-trivial potentials and kinetic terms, etc.. All these 
C/3 models have to pass the observational constraints imposed by current data on the cosmological background as well 
, as on the power spectrum of the perturbations. In many cases the lack of other observables renders many models 

observationally indistinguishable. However some models have other observables like non-Gaussianity (i.e. non-linearity 
y—i ^ of the perturbations) and these can be and are being used to further constrain and distinguish models. 
^> . Non-Gaussianity of the primordial curvature perturbation has not been observed, however the CMB data still 
On ' allows a small non-Gaussianity (~ 0.1%) and there have been some claims of detections or hints of its presence (sl-flH. 
Despite this small value, efforts are underway (e.g. ESA's Planck satellite [l2|) to measure the CMB anisotropics to 
a precision with which this small amount would be detected to many standard deviations. 

It also turns out that the simplest models of inflation, i.e. the ones driven by a single scalar fleld with standard 
J^j ■ kinetic term, with a potential satisfying the slow-roll conditions and if one assumes the standard initial conditions 
' for the quantum vacuum fluctuations, predict a small and unobservable level of primordial non- linearity |13H18|. 
I Therefore only models that can produce large non-Gaussianity can be constrained if no detection is made. One of 
• • ' such models put forward by Starobinsky in 1992 [l^ will be the focus of this work. Many models which break at 
^ , least one of the above conditions have been proposed and studied recently. An incomplete list is: for models that use 
k> • several flelds [20-51], for models with non-canonical kinetic terms (53 - [63| . models with temporary violations of the 
j_] , slow-roll conditions or small departures of the initial vacuum state from the standard Bunch-Davies vacuum |63 - [73j . 
I For reviews on the theory of non-Gaussian perturbations see ^73-[78| . 

The initial motivation of the Starobinsky model [19] that breaks temporarily the slow-roll approximation was to 
explain an observation that the correlation function of galaxies seemed to require more power on large scales than 
at smaller length scales within the paradigm of a Einstein-de-Sitter universe (a spatially flat universe with cold 
dark matter particles). It was also used to explain the possible nonconservation (enhancement) of the curvature 
perturbation on superhorizon scales due to not-fast-enough decay of the decaying mode [t^. More recently similar 
models [l^ [sol - issj were revived and proposed because it was realized that they can provide better flts to the CMB 
anisotropics power spectrum than the primordial power-law model [86l - [9l| . Another reason for this interest is that it 
was realized that such models have very distinctive bispectrum signatures like for instance non-linearity parameters 
that are strongly scale-dependent. A particularly appealing feature of the Starobinsky model, a vacuum energy 
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dominated potential added with a linear term that changes slope abruptly at a point, is that despite breaking a 
slow-roll condition it admits an analytical treatment of the background and linear perturbations equations as well as 
the non-Gaussianity as was shown in [tO, [Q^I and in the present paper. 

The main goal of this work is to compute analytically the leading-order bispectrum for the Starobinsky model. The 
bispectrum on large scales was first computed in ^7Q] using a next-to-leading order gradient expansion formulation 
and more recently on all scales in the equilateral limit in [92] , we compare the results of our calculation with these 
works. 

This paper is divided in four sections and an appendix. In Section |TT1 we describe the model and the analytical 
solutions for the background and linear perturbations are briefly reviewed. In Section IIIIl we compute the leading- 
order contribution to the bispectrum in the equilateral limit. In the second subsection we shall discuss and compare 
our results with previous attempts in the literature. In the last subsection we elaborate on the size and shape of the 
analogue of the non-linearity parameter in this model. Section lTVl is devoted to conclusions. In Appendix[X]we present 
the lengthy result for the dominant contribution to the bispectrum that is valid for any triangular configuration (not 
only in the equilateral limit) of the three momenta on which the bispectrum depends. 

We work in units where the reduced Planck constant h and the speed of light c are h = c = 1 and the reduced 
Planck mass is Mpi = (SttG)"^/^, where G is Newton's gravitational constant. 



II. THE MODEL 



In this section wc introduce the model under consideration and review the analytical solutions for the background 
and the perturbation equations. All the results of this section can be found in Refs. 19, 92], here we only briefly 
review the main results that will be needed in the computation of the bispectrum of the next section. We use a 
somewhat different notation at points. 

The Starobinsky model consists of a canonical scalar field minimally coupled to Einstein gravity with a potential 
given by 

where (j) is the scalar field value, Vq, A+ and A„ are model parameters that are assumed to be positive. The 
sharp change in the slope occurs at = 0o and (po is another free parameter of the model. We also assume that 
\y — Vo\/Vo ^ 1, which means that the potential is vacuum energy dominated. 

In the present case we are interested in a Friedmann-Lemaitre- Robertson- Walker flat, homogeneous and isotropic 
background spacetime given by 

ds^ = -dt^ + a^{t)5,jdx'dx\ (2) 

where t denotes cosmic time and a{t) is the cosmic scale factor. The relevant equations of motion are the Friedmann 
and Klein-Gordon equations given respectively by 

where the Hubble rate is defined as usual H = a/a and dot denotes the derivative with respect to t. The slow-roll 
parameters are defined by 

(4) 




and will be assumed to be small before the transition. After the transition e will continue to be small but rj will 
become large for some time before the slow-roll regime is again recovered. 

With the above mentioned assumption of a vacuum dominated potential the dynamics for the scale factor becomes 
particularly simple and is described by 



(5) 
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which has a solution a(t) = exp(_ffoi), where Hq = ■\/Vo/(3M|>;) and we chose a{t = 0) = 1. In conformal time, 
denoted by r, the previous solution is a(r) — —l/{HnT ) an d one has —Hot = \n{—HaT). Before the transition the 
Klein-Gordon equation can also be easily solved to find |19l . lo^ 

0+ ~ 0. + + ^' In(-ffoT), -P- ~ -aH + , (6) 

where 4>i is the initial value of the field. The subscript + denotes quantities before the transition and the subscript — 
denotes quantities after the transition. 

After the transition, the Klein-Gordon can also be solved despite the fact that the slow-roll parameter rj- becomes 
temporarily large (however e_ is small), one finds [l^l 

(7) 

where AA ~ A^ — A^ and tq is the transition conformal time. The slow-roll parameter e before and after the 
transition can be written respectively as 

where = —{AA/ A_)kQ and kg is the transition scale given by /cp — o-{tq)Ho = — I/tq. It is worth noting that 
before the transition ry+ ss 4e+ and at late times (much after the transition) we have again rj- w 4e_. For a transient 
period after tq, rj- can become large and it is well approximated by 

~ ~ (9) 




1 — p'^T 



,3^3 ■ 



In Ref. 92| it has been shown that the previous expressions for the slow-roll parameters agree very well with the exact 
numerical solutions. For the plots of the next section we choose the same parameters as the main choice of Ref. [9^ . 
i.e. Vo = 2.37 x lO-i^Mj,;, A+ = 3.35 x 10-"M|,,, A_ = 7.26 x IQ-i^Ml,;, this gives AA/A_ w -3.61. This choice 
satisfies the constraint on the normalization of the primordial power spectrum imposed by COBE observations. 
We now briefly describe the analytical approximation for the mode functions, all the details of the derivation can 



be found in H, . 



We work in the comoving gauge, so that the background value of 0, i.e. (j)o, determines when the transition happens 
even including perturbations. In this gauge the three-dimensional metric is perturbed as 

= a^e'^S.j, (10) 

where TZ denotes the comoving curvature perturbation and tensor perturbations were neglected because they do not 
contribute for the tree-level scalar bispectrum. 

The Mukhanov-Sasaki equation in Fourier space is 

n'l + 2-7^'fe + = 0, (11) 

z 

where z = a\f2k, prime denotes derivative with respect to the conformal time t and the Fourier transform is defined 
by 

7^k = 7^(r,k)= J d^xn{T,x)e-'^ ''. (12) 

The standard quantization procedure in quantum field theory is to promote TZ to an operator that is expanded in 
terms of creation and annihilation operators and mode functions as 

TZ{t, k) = TZ{t, k)a(k) -f- 7^*(T, -k)a^(-k), (13) 

where the annihilation operator a and the creation operator satisfy the usual commutation relation [a(k), a^^(k')] = 
(27r)'^i5^'^' (k — k'). The power spectrum of the curvature perturbation is given by 

27rP^\k) ^ V2k^\nk{Tf)l (14) 
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where t/ is the end of inflation time. 

Before the transition, the slow-roU approximation is vahd and the mode function solution of Eg . (1111) t hat reduces 
to the usual Minkowski result on small scales is the well known leading order slow-roll expression [l^ [921 ; 



7^+(r,k) 



and one can take the derivative to find 



n(r,k) 



3gg 1 

A, 



1 



(fcr — i) e 



A, 



(15) 



(16) 



This last equation is valid at leading order in slow-roll and will be used in the computation of the dominant contribution 
to the bispectrum. After the transition an analytical approximation can also be found, see [l^, for all the details. 
Here we just present the final result which is given by 



7e_(T,k) = 



iHoak 



where the Bogoliubov coefficients are 



(1 + ikr) e 



(1 - ikr) e 



(17) 



ttfc = 1 



iil\A ko 



1 + 

k^ 



/3fc = - 



MAA ko 
2A+ T 



k 



(18) 



The most important feature of the above solution is the presence of the negative frequency mode which introduces 
oscillations in the power spectrum. For the derivative one finds 



7e'_(T,k) = 



iHaak 



-He_{l + ikT) 



(1 -I- ikr) + k^T 



2Mpiy/We 



-Ue- (1 - ikr) - 



Hf]- 



(1 - ikr) + k^ 



-ikr 



AkT 



(19) 



where % ~ aH is the conformal Hubble parameter. This expression will also be used in the computation of the 
dominant bispectrum contribution in the next section. Because we are only interested in the dominant contribution 
we neglect the first term inside the square brackets because it is proportional to e which is always small in the present 
model. 

III. THE EQUILATERAL BISPECTRUM 

In this section we compute the leading-order dominant contribution to the bispectrum in the equilateral limit. 
After that we compare our result with previous results in the literature. Finally we discuss the shape and size of the 
analogue of the non-linearity parameter in this model. We leave for Appendix |X] the lengthy result of the dominant 
contribution to the bispectrum that is valid for any triangular configuration of the momenta. 

A. Calculation using the two vertices action 

In order to compute the tree- level scalar bispectrurn, we use the in-in formalism (93l . |93 | . For this we need to find 
the cubic-order interaction Hamiltonian, see e.g. Ref. M for a review about this procedure. 

The third order action after ignoring total derivative terms can be found for instance in fl^, [sl, [9^ [9^ . It reads 
(here we follow the notation of 96]) 



Mpi I dtd^x 



a^e di] 



2a 



\ SL 




) sn 


1- 



(20) 
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where we have defined 



5L 
5n 



V dt 

-{dTZ){dTZ) + d-^id^djid^-RdjTZ))] 



1 



(21) 



[[dn){dx) - d-\d,d,{dmjx))\ , (22) 



and d^^ denotes the inverse Laplacian. 

In the present model, the slow-roll parameter e is taken to be always small after and before the transition, this 
means inflation never stops. However both rj and rj' can become large temporarily after the transition. Therefore the 
dominant contribution to the bispectrum is expected to come from the following terms in the action, 



S. D M 



dtd'^x 



2 dt 



sn 



+ Ml, / dtd^x 



dt 



2 



(23) 



where the last term is a total time derivative but it is retained because it is proportional to jLThis is one of the 
many terms that appear when simplifying the action to the form (j20p via integrations by parts [qtI [98| , which can be 
important and cannot a priori be neglected [9^ . After integration by parts the cubic action can be simplified to 



S3 3 Mpi 



dtd^3 



from where one easily finds the cubic-order interaction Hamiltonian as 



Hint (t) 



Mil 



The tree-level three-point correlation function (or bispectrum) at the time Tg after horizon exit is 



dra(0| [n{Te, ki)7e(Te, k2)7t(re, ka), iJ„t (r 



(^]|7^(re,kl)7^(Te,k2)7^(re,k3)|^^) = - 

More explicitly, the bispectrum is 

(n|7^(0,kl)7^(0,k2)7^(0,k3)|^]) « {2TTf6^^\K)2M 



(24) 



(25) 



(26) 



7^(o,kl)7^(o,k2)7^(o,k3 



(27^'*(T,k2)7^'*(T,k3) - kln*{T,k2)n*{T,k3)) 



dTT]ea^n*iT,ki) X 

two perms., (27) 



where K = ki + k2 + k3 and "two perms." denotes two additional permutations of the term displayed, and in the rest 
of this work we set Tg « 0. 

Because the integrand of the previous expression is different before and after the transition, here we compute the 
contribution before the transition, i.e. the integral from —00 to tq = — 1/fco, separately from the contribution after 
the transition, i.e. the integral from tq to zero. Obviously the final answer for the bispectrum is the sum of these two 
contributions. As we show below, for the parameter choice described in the previous section, the contribution after 
the transition is much larger than the contribution before the transition, so the later one can be ignored. 

The contribution before the transition (i.e. the integral is performed from —00 to tq only) is 



(^^|7^(re, kl)7^(r„ k2)7^(r„ k3)|r!)_ 



-(27r)3,5(3) (K) 



1 



32M 



2z k- 



21,2 



2 '^3 



e [ 1 + — + I — 



+kle ( fc, 



fcifc2^3 .kik2k3 .kik2 + k2k3 + kik3 



ktko 



+ i- 



two perms. 



(28) 



where kt = ki + k2 + k^. In the equilateral limit and after using the explicit expressions for the slow- roll parameters 
it simplifies to 



97 ff ® A 1 

(^]|7^(o,k)7^(o,k)7^(o,k)|^^)+ = -{2T:f5^^^ (k) 4 , .° t — 



4 M2;A3_ fc6 



{ak-PkTe~^'T^^ [2i 



k ko 
ko k 



(29) 
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In the large scales limit, i.e. ^ fco, it further simplifies to 

{n\n{o, k)n{o, k)n{o, k)\n)+ 



97 FfS 

(27r)3j(3)(K^^^ ^0 



4 Mf,iAl fc6 



O 



f 1 



(30) 



The previous bispectrum is at leading order proportional to This is the scaling expected for a scale invariant 

model. 

On small scales, i.e. when fco ^ k, we get 



— sin ( ^ I + ( 2 - 



27 H^A+ 1 

217/ 



cos 



9AA /^M^ro/^^o 



(31) 



where the first term in square brackets is the leading one and the others are the sub-leading corrections, one of 
which oscillates with a different "angular frequency". As one can see, in this limit the bispectrum scales as 
This represents a linear growth envelope with respect to the bispectrum of a scale-invariant model. The "angular 
frequency" of the leading-order term is set by the transition scale and is S/feg. 

The calculation of the contribution after the transition (i.e. the integral from tq to Te) is more involved but can be 
done analytically. The full result can be found in Appendix [Xj In the equilateral limit it simplifies to 



(^!|7^(o,k)7^(o,k)7^(o,k)|^^)_ = -{2nf5^^^K'^''^^'' ^ 



xSR 



2 At fc9 



where the functions Fi and F2 are defined as (for Te 

2 



0) 



Fi 



p^i 



- 1 -e 



-3! 




k^ 



;„5 



(32) 



(33) 



'3i+{ — 



9p^ 



6k''^p^ k^p^ 



1,5 



where p^ — —{AA/Aj — /cq. Equation ([5^ is one of the main results of this paper. 
On large scales, i.e. A; ^ fco, it simplifies to 



(fi|7t(0, k)7t(0, k)7^(0, k)|f7)_ 



-(271)3(5(3) (K).^'^o'^^ ^ 



20 A^Ai fc6 



41 4A+ 
70 ~ 5AZ 



k 
fco 
31/12 



jL6 



(34) 



-38851 4246Ai 



113400 14175A_ 



AAl 
25A^ 



ky f 225409 5701/1+ ^ 

fci^y 1^3492720 3969MI ~ 2835^2; 



(35) 



where the first term in square brackets is the leading one and we have included several other sub-leading terms 
because they will be used to compare our result with previously known results available in the literature. On large 
scales the quantity fc^ x bispectrum is proportional to tt so its amplitude decreases for larger and larger scales 
until this contribution becomes smaller than the contributions coming from the slow-roll neglected terms in the 
previous calculation. These contributions are expected to be small and perhaps unobservable. Some (but not all) of 
these contributions were computed in (9^ and they turn out to be of the scale-invariant kind (that is the quantity 
fc^ X bispectrum is independent of fc). 
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On small scales, i.e. when fco <iC k, we get 

06 /t12 A /I 1 

(^^|7^(o,k)7^(o,k)7^(o,k)|o)_ = (2^)35(3) (k) 



4 k^ 









X 













' 2 - -— cos — - -— cos — + C 



A+ 2A+ J \ko J 2A+ \ko J \ k 



(36) 



In this limit, the leading order term (first term inside the square brackets) is a linear growth envelope with oscillations 
of an "angular frequency" 3/fco. This linear growth may potentially have important implications regarding the 
detectability of this signal. This result is similar to the one found for instance in Ref. [o^. Also for the parameters 
values choice of the previous section and used in [o^ one can show that the contribution (IM]) to the integral dominates 
over Eq. (j3D). In fact the ratio of these amplitudes is proportional to AA/{A+e+) which is always large for the case 
of interest. Similar result holds on large scales, where Eq. (l35t dominates over Eq. (|30|) (this is true for scales larger 
than the transition but not much larger, because Eq. ([55t decays quickly as the length scale increases). 

Finally, it is instructive to estimate the range of wavenumbers Afc over which large deviations from Gaussianity 
maybe be expected in this model. If the sharp transition of the Starobinsky model ([T]) is replaced with a smooth 
transition with a width in field space then the number of e-foldings that takes to cross the transition is AA^ ^ 
i7o/0A0 ~ -Vb/(M|,;A_)A0. Because Afc/fco l/|AiV| one finds Afc/fco (M|,;A_)/(A(?!)Vo). In the present 
model A0 is zero which implies that the range of scales affected is Afc oo. In a more realistic model, we expect 
that any transition is eventually smooth in a sufficient small scale, this transition width will introduce a cut-off for 
the growth of the bispectrum as described by the previous equation. For scales smaller than this cut-off scale the 
previous equation is not valid, instead we expect the amplitude of the deviations from Gaussianity to go quickly to 
zero. This expectation was recently confirmed in a different but somehow related model (the potential has a step-like 
feature) 



B. Comparison with previous works 

As mentioned in the introduction the bispectrum in this model was first computed in fl^. There, the spatial 
gradient expansion approach was used at next-to-leading order. Because of this their result is expected to be correct 
only on large scales and also in the limit A+ ^ In the equilateral limit, the first equation in Eq. (6.20) of 7^ is 



(^7|7^(o,k)7^(o,k)7^(o,k)|^]) = -(271)^5(3) (k; 



3^ H^^AA 1 
10 A^AX fc6 



kV 4/A:\V. A+\ 4 /fc\V. .A+ Al\ I f k V f , .A+ A^ 



koj ' 5\koJ A^)^25[ko) V ^A^+A^j + gUoj V '^'A^ ' A^ 



(37) 



Using the fact that A^ ^ A^ the above equation is, up to a factor of 2, equal to Eq. (|35|) at leading order but it differs 
in the term proportional to {k/ko)^. This extra factor of 2 in Eq. ([57| is due to a double counting mistake of the 
permutations in ,70J^. We have not found a reason for the discrepancy in the term proportional to (fc/fco)*. Finally 
one can compare the result of [tO] for a general triangular configuration with Eq. (jAl[) of Appendix [A] Correcting 
the factor of 2 difference in the overall amplitude due to the double counting mistake, we again find that their result 
agrees with ours up to a discrepancy in the sub-leading term of the order of fc^/Zcp. 

The second attempt [ol] to compute the bispectrum (equilateral configuration only) for this model used the in-in 
formalism just like in the present paper and even computed analytically some of the sub-leading order corrections 
(although not all of them). They found that the dominant bispectrum contribution on both large and small scales 

goes as 1/k^ and for small scales the amplitude oscillates with cos (^f^)- discussed in the previous subsection we 

found a different result. On large scale the bispectrum scales as l/(fcgfc'') while on small scales it possesses a linear 
growth envelope as l/(fcoA;^) sin(3fc/fco). We discuss the origin of this difference below, but in short it appears that 
Ref. missed a term in the calculation that is proportional to the Dirac delta function. 



^ We thank Yuichi Takamizu for discussions on this point. 
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In Ref. [9^ an alternative form of the action^ given in Eq. ([23)) is used. Keeping the term proportional to 577' only 
one finds 

(^^|7^(re,kl)7^(Te,k2)7^(T„k3)|^^) = (27r)3,5(3) (K)^Af|,^7^(Te,kl)7^(Te,k2)7^(Te,k3) 

draper]' (7^*(r, kl)7^*(T, k2)7^*'(T, kg) + two perms.) +c.c., (38) 



where "c.c." denotes the complex conjugate. The action (l23l) contains two further terms, one proportional to the 
first order equations of motion, but this term gives an exactly zero contribution when evaluated on-shell. The second 
term is a boundary term. Which should be evaluated at the past infinity boundary, this gives zero once we rotate the 
contour into the imaginary plane to select the free vacuum. And it should also be evaluated at a later time (a time 
when all scales of interest are well outside the horizon). However at such a late time after the transition, slow- roll will 
have been restored and again the contribution from this term is sub-dominant. Therefore one concludes that Eq. 
contains the dominant contribution. 

Using the Klein-Gordon equation, 77 can be written as 



V 



2 cfV 



2 A, 



H |0'(ro)| 



■S{t - To) + 



(39) 



where • • • denotes other terms without Dirac delta functions, the last equality is valid for the Starobinsky model and 
we have used S{(f>o — (j)) = S{t — tq) /\<j>' {to)\. If we keep only this term in Eq. ([55]) the integral simplifies to 

(^^|7^(Te,kl)7^(re,k2)7^(re,k3)|^!) = (27r)3^(3) (j^^ 



4M2je(To)(A+-A_)a3(To) 



l'^'(To)|i?(To) 

x5 [7^(re, ki)7e(Te, k2)7e(re, k3)7^*(ro, kl)7^*(To, k2)7^*'(ro, kg)] 

-I- two permutations. (40) 

This contribution dominates over the dominant contribution calculated in (92| on small scales and is of the same order 
of magnitude for large scales. 

For the equilateral configuration, i.e. ki — k2 — — k, the bispectrum (|40|) with Tg = simplifies to 



(^]|7^(0,k)7^(0,k)7^(0,k)|^7) = {2T:f6^^^^ (K) — 



AIA^_ k^ka 



-■ii-r- 

e '-0 



In the large scale limit, i.e. k ko, this reads 

(^7|7^(0,k)7^(0,k)7^(0,k)|^7) = (27r)^,5(^) (K 
On small scales, i.e. ko ^ k one finds 

(^^|7^(o,k)7^(o,k)7^(o,k)|^]) = (271)^,5(3) (k) 



AAH, 



12 r 







^3fc\ 


X 










Kko) 



4 J Al fc6 



3^^A_H^ 

" 2A+ 



UA^ - 2A^ 

ioaZ 



O 



4n 



(41) 



(42) 



9AA 
2A~ 



(43) 



If one calculates Q /k^ ~ fc^ x bispectrum (the precise definition of G is given in the next subsection) then on small 

scales one finds that G/k^ ^ sin which is the linear growth found in and in the previous subsection. On 

large scales G /k^ ^ constant. In order to obtain the final answer for the bispectrum one needs to add to the above 
the contribution from all the other terms that was first calculated in Ref. 193. 



One is free to use this form of the action or the form used in the previous subsection or possibly others, the end result is the same as 
we shall show below and the two procedures are equivalent. 
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Using their Eqs. (100)-(104) one has that the contribution is 

(^7|7^(o,k)7^(o,k)7^(o,k)|^2) 



(2^)3^(3) (K) 



^i(fc) sin — 
ko 



4 fc6 
k 



A2{k) cos — 



3fc 

yl3(fc) sin ( — ) + A4ik) cos 



where 



A2{k) = 

Asik) = 
Aiik) = 



3A^_AA 
AA% 

3A^_AA 



12A\_ 



Al 
12A5 



k^ 



9A- 

2 r 



1 + +A 



9A_ 



k^ k^ 



-A, 



27(AA)2 1 - - 27AA (5A_ - 7A+) - (9A_ - 11^+)^ 



- 27^2 
+Al ( 81 



fc2 

^0 



-2 I +18A_A4 



/fc2 fc4 fc6 ^^8 

153^-9^-93^+4^ 

fi,Q h,Q /iQ /VQ 

On large scales, i.e. A: <C fco, one finds that the bispectrum (|44|) reads 

56 A 4 ffl2 



fc2 



A:2 fc4 



(^]|7^(o,k)7^(o,k)7^(o,k)|^^) = {2n)H^'^ (k) 



1 + 



17A_ - 
10v4_ 



On small scales, i.e. ko <^ k one finds 



(^l|7^(0,k)7^(0,k)7^(0,k)|^]) = (27r)3(5(3) (K) 



36 AA r 



cos 



3fc 



3fc' 
^0 . 



,(44) 

(45) 
(46) 

(47) 



(48) 



(49) 



(50) 



Then it can be seen that the sum of Eq. ([^^ with the Dirac delta function contribution given by Eq. (PT|) exactly 
agrees the result in Eq. ([5^ given in the previous subsection. This explains the discrepancy between our result and 
the result of [o?] (it also explains the previously existing discrepancy between the large scales results of fo?] and ) . 
The point is that the Dirac delta function term in rj of Eq. (jJH) cannot be neglected because on large scales it gives 
a contribution comparable to Eq. (HHl) and on small scales it actually dominates over Eq. ([5(7)). The plots of these 
different contributions for the final bispectrum are given in Fig. [TJ with the choice of model parameters described in 
Section HH 



C. The non-linearity function Q{k)/k^ 

In this subsection, in order to gain a feeling for the shape and size of the bispectrum in the equilateral limit for the 
Starobinsky model, we plot the previous results for the choice of parameters as in Section im 

To show the shape of the bispectrum when the three wavenumbers are of comparable size one commonly plots a 
quantity defined by 

- ("l«(0.M«(0,.,)m.3)|O,, ,51, 



If this quantity is of order unity there may be some hope that it will be measured in the near future for example with 
ESA's Planck satellite. Because in the present model the amplitudes of the power spectrum before and after k^ can 
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|BispectruiTi| 




0.01 0.1 1 10 100 

k 



FIG. 1: Plot of the absolute value of the equilateral bispectra for fco = 1 as a function of k. The red long-dashed line is the 
plot of Eq. (|44[) . the green dotted line in the plot of the Dirac delta function contribution to the bispectrum as given by Eq. (|4ip 
and the solid black line is the plot of Eq. (|32|) . One can see that for large scales the contributions from Eq. (|44|) and Eq. (|4ip 
are comparable in size and because they have opposite signs they partially cancel to give the solid black line. On small scales 
the amplitude of the long-dashed red line is negligible compared with the dotted green line. For this scales the dotted green 
line and the solid black line are virtually identical. The model parameters have the values described in Sec. |lll 



be quite different it makes more sense to define two functions 5< and t/> to study the bispectrum for scales larger 
and smaller than the transition scale fcp respectively, as 

'-^^ - (o|«(o,.OK(o,.,)K(o,.3,|a)., ,53, 

where P< and P> denote the asymptotic values of the power spectrum from large and small scales respectively, and 
are given by 

With this last definition we make sure that the oscillations seen in the quantity are truly bispectrum oscillations 
and are not due to the oscillations already present in the small scales power spectrum. 

As one can see from Fig.[5J on large scales and with our choice of parameters of Sec. [IT] the amplitude of G</k^ is 
small and decays quickly towards larger scales. At some point this amplitude becomes so small that it is subdominant 
with respect to the slow- roll correction ignored in this work. Some of these slow-roll corrections were computed in fo^ l . 
On the other hand. Fig. [3] shows a very interesting behaviour for the non-linearity parameter G>/k^- One finds rapid 
oscillations with a growing amplitude towards smaller scales. For the model parameter of Sec. |lT]the amplitude can 
easily become much larger than one. 



IV. CONCLUSION 



In this work, we revisited the computation of the tree-level three-point function (or bispectrum) in the Starobinsky 
model of infiation. In this model, a vacuum energy dominated potential is supplemented with a linear term in which 
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10-5 r 

10-': 

0.01 0.05 0.1 0.5 1 

k 

FIG. 2: Plot of Q^/k^ as a function of k for ko — 1 and the model parameters of Sec. [Ill On large scales one finds that the 
amplitude of Q</k'^ is much smaller than one and is decreasing for larger scales. This is an expected result as for these scales 
the transition happened when they were super-horizon, because TZ is non-linearly conserved the deviation from Gaussianity 
should be small as these perturbations were nearly Gaussian at horizon crossing. For very large scales, eventually the magnitude 
of this contribution becomes smaller than the always present slow-roll corrections that we neglected in this work. When this 
happens these sub-leading corrections become the dominant part. However they are slow-roll suppressed and so less interesting 
observationally. 




FIG. 3: Plot Q->/k^ as a function of k for ko — 1 and the model parameters of Sec. |lTl On small scales the amplitude of 
the analogue of the non-linearity parameter, i.e. t/>/fc'^, oscillates with an "angular frequency" given by 3/fco, where fco is the 
transition scale and the envelope of the oscillations grows linearly with fc/fco. One easily sees that the amplitude of Q^jk? can 
reach values much larger than one and be of potential interest observationally. 

the slope changes abruptly at a point. It is well known that under the assumption of vacuum energy domination this 
model admits an analytical solution not only for the background equations of motion but also for the mode functions. 
Thus despite the fact that the slow-roll approximation is broken for a period of time just after the transition happens 
this model still admits an accurate analytical approximation for the mode functions of linear perturbation theory. 

Using this analytical approximation for the mode function we computed analytically the dominant contribution to 
the bispectrum. The main result is Eq. (j32p . In the equilateral limit, we obtained that on scales larger than the 
scale of the transition fco the dominant contribution to the bispectrum scales as l/(fcQfc^), that is, the analogue of the 
non-linearity parameter, Q/k^, scales as Q/k^ ^ (fc/fco)^. This is in agreement with the leading-order result of [t^ 
(after correcting its factor two error) but it differ from the large scale result of [92,] . This last result was found to scale 
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as Q /k'^ ~ constant. We explained the origin of this difference with the fact that Ref. 'o^] missed the contribution 
coming from a term in the cubic-order action that contains the Dirac delta function in the time derivative of the 
coupling constant, efj. This missing contribution has the opposite sign to the dominant term they calculated and it 
partially cancels it out and one is left with the result we found, that is Q/k^ ~ {k/ko)'^. For small scales, we found 
some interesting behaviour. The parameter G/k^ scales as G/k^ ~ k / ko sm{3k / ko) . The "angular frequency" of the 
small scale oscillations is set by the transition scale and is given by 3/ The envelope function of these oscillations 
in a linear growth toward small scales. A similar result was found in '9&\ for a somehow related model. This linear 
growth represents a strong scale dependence of the non-linearity parameter in this model and the k/ko enhancement 
factor can have important consequences regarding observations. It would be interesting to consider observational 
constraints on this kind of oscillating and growing amplitude G/k^ models, however that study is beyond the scope 



of the present work. Once again this result differs from the result found in Ref. [92| . There, it was found that G/k^ 
scales as C//fc^ ^ cos(3fc/fco), with no enhancement factor. The reason for the difference with respect to our result is 
again the fact that the contribution from the Dirac delta function was missed. Finally, we briefly showed that the 
amplitude of G/k^ can reach values much larger than one on small scales. The large enhancement factor k/ko is key 
in reaching that conclusion. 

For future work, we leave the computation of the trispectrum or four-point function. 
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Appendix A: The bispectrum for any configuration of the triangle 



In this Appendix we will present the result for the bispectrum that is valid for any allowed configuration of the 
three momentum vectors. 

As discussed in the main text, the integral in Eq. (I27p can be split in a contribution before the transition and a 
contribution after the transition. For a general triangle configuration, the contribution before the transition is given 
by Eq. ([SS)) and it is sub-leading with respect to the contribution after the transition. Therefore for the scales of 
interest the bispectrum is well approximated by the contribution after the transition only. It turns out that the 
integral after the transition can also be evaluated analytically. After a lengthy but straightforward calculation, the 
result for the bispectrum is 



(^^|7^(Te,k)7^(re,k)7^(re,k)|^])_ = -(2^)3,5(3) (k) 



x3? 



("fei - Pkt) (afc2 - Pk2) ("fca - A-a) {^*ki^*k2^*k3T{ki,k2,k:i) - l3l^(3l^f3l^T{~ki, -k2, -fca) 
-(^l^otl^PlT(ki,k2, -ks) - al^al^l3^^T{ki, -fcz, fcg) - al^al^(3l^T{-ki,k2, fcg) 
+<Pkl3knki,^k2, -fcs) + alPlPin-ki, k2, -k,) + all3l/3lT{~k,,~-k2, k,) 



(Al) 



where T{ki, k2, k^) is defined as 

„3 



T{ki,k2,k3) = 



Jrkt 



(1 - p3r3)3 



- 3 + kik2k-iktp^T'^ + ip^T^ (9fcifc2fc3 + ki{kl + kl) + k2{kl + kj) + kaikj + kj)) 
-p^T^ (kj + kl + kl + 9{kik2 + A;2fc3 + ^3^1)) - (fcifc2fc3^t + Qip^h) 

{9p^ - ikt{kik2 + k2k3 + ky,ki)) + {kl + kl + kl + 3(fcifc2 + fc2A:3 + ^3^:1)) + iikti 



(A2) 
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For Te = 0, the function T{ki, ^2, fca) is related to the functions Fi and F2 defined in the main text as Fi = T{k, k, k)/6 
and F2 = —T{k, k, —k)/2 and in the equilateral configuration Eq. (|Aip reduces to Eq. ([5^ . 

Finally it is worth noting that in the squeezed limit of the triangle formed by the three momenta, i.e. when 
^2 ^ fca ^ fc and the small momentum ki is ki <^ k, one can find on large scales, i.e. when fco S> fc 3> fci, 

(o|7^(o,k)7^(o,k)7^(o,k)|^7)_ « -{2nfs('^ (k) ^^l^^ (^Aj , (ab) 

and on small scales, i.e. when ko ^ ki ^ k, one finds 

(f]|7l(0,k)7t(0,k)7t(0,k)|f^)_ « (2.)3<5(3) (K) f sin g)) , (A4) 

where the enhancement factor fci/fco in the previous equation is fci/fco 3> 1. 
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